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Abstract 

If the scalar tau f is the next lightest supersymmetric particle and decays into a gravitino(G) 
being the lightest supersymmetric particle, it will have generally a very long lifetime. In this paper, 
we investigate the possibility to study the decay of such a long lived scalar tau at the LHC. If we can 
add to the present LHC experiments additional detectors which are able to stop the stau particles 
and measure the produced decay products, the decay characteristics can be studied precisely at 
the LHC. We identify a maximum "stopper detector" that could be added in the CMS cavern, and 
estimate the sensitivity to the lifetime of the stau and to the mass of gravitino with this detector. 
The decay of the scalar tau may be significantly modified if the decay channel to the axino a is 
open. We study the possibility to distinguish such decays from decays into gravitinos by measuring 
the process f — > cl(G)tj using the stopper detector. 

PACS numbers: 
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I. INTRODUCTION 



The minimal supersymmetric standard model (MSSM) is one of the most important can- 
didates for the physics beyond the Standard Model. The strongly interacting superpartners, 
gluinos and squarks with mass lighter than 2.5 TeV can be discovered at the Large Hadron 
Collider at CERN (LHC) fl. The physics run of LHC will start in 2008. 

Among the supersymmetric (SUSY) particles, the lightest SUSY particle (LSP) plays a 
key role. In cosmology, it is a candidate for the cold dark matter in the universe. At the 
LHC, the signals of the supersymmetric particles depends on the nature of the LSP. It may 
be the lightest neutralino Xi> which escapes from detection and leads to missing E? in the 
event. Another possibility for the LSP is the gravitino G, the superpartner of the graviton. 
The gravitino coupling to other particles in the MSSM sector is extremely small, namely 
proportional to 1/M pl . The next lightest SUSY particle (NLSP), whose decay products 
necessarily include the gravitino, is therefore long-lived. 

We discuss the scenario where the LSP is the gravitino and the long-lived NLSP carries 
charge. A natural candidate for a charged NLSP (CNLSP) is the lightest scalar tau, f 1; 
which can be significantly lighter than the other sleptons due to the left-right mixing induced 
by the off diagonal matrix element m\ R ~ m T /xtan/3. The charged particle leaves tracks 
in the central detectors (ATLAS and CMS), giving additional information for the SUSY 
particle reconstructions. If the NLSP decays in the main detectors, a displaced vertex may 
be observed as well. 

The expected lifetime of the CNLSP f is unconstrained, because it is proportional to 
(tUq) 2 the yet unknown gravitino mass squared. On the other hand, the gravitino mass is 
proportional to the total SUSY breaking scale in the hidden sector, therefore the determi- 
nation of the lifetime is an important physics goal. The lifetime measurement gives direct 
information on the hidden sector. 

A particle decays efficiently in the main detectors (CMS and ATLAS) if the decay length 
is sufficiently short, cr <C (10m)x A" pro( iuced, where Adduced is the number of produced SUSY 
particles. On the other hand, for typical SUSY production cross sections, a direct observation 
of the decay is very difficult for tcnlsp > 0.01 sec. However, it has been pointed out that 
the CNLSP stopped by the ionization .os 8 in the materia, may be studied in detai. fl g 
An idea for a stopper based on a water tank is presented in Ref.[3|, where the water can be 
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transported away from the detector site, for concentration and further study. In Ref. a 
detector consisting of a tracker and heavy stopping material is proposed, which can measure 
the arrival time and the location where the CNLSPs are stopped, in addition to the energy of 
the decay products. Another possibility, which requires minimal experimental modification, 
is to study the CNLSP which are stopped in the main detector or surrounding rock [J]. 

As pointed out in Ref. [5|, the study of the CNLSP decays can probe the underlying 
supergravity in nature. With the gravitino mass inferred from the kinematics, the additional 
CNLSP lifetime measurement will test an unequivocal prediction of supergravity. Moreover, 
the study of a rare 3-body decay f — > t^G can reveal the peculiar couplings of the gravitino 
and the gravitino spin 3/2. 

In this paper, we consider the physics of the CNLSP f decays that can be done with the 
heavy material stopper-detector 2|, because only a detector of this type can cover a wide 
range of lifetime O(10nsec) < t CNL sp < O(10 years). We find that the mass of the gravitino 
can be measured if it is sufficiently heavy (roughly ttlq > 0.2rrif). In that case, one can 
check if the lifetime is consistent with the supergravity interpretation. 

As we shall see, the LSP mass resolution is however poor if itlq < 0.2mf . In that case, 
it is very hard to prove that the decay f — > r plus an invisible particle X is indeed caused 
by the supergravity interaction involving G. We should note that, because of the extremely 
weak coupling of the gravitino, if there is any unknown supersymmetric particle X which 
couples rather weakly to the MSSM particles, the lightest MSSM particle may decay into 
X instead of G even if the G is the LSP. For example, the superpartner of the axion, the 
axino a, can be such a particle. However, we found that it may be possible to discriminate 
the case of X = G and the case of X = a, by investigating the three body decay f — > t^X 
as suggested in [5|, |6j, with enough statistics. 

This paper is organized as follows. In section [Til we discuss a possible detector setup. We 
found an 0(1) kton detector (up to 8 kton) may be placed next to the CMS detector without 
serious modification of the CMS experiment itself, but with non-negligible modifications to 
the CMS cavern side walls. In section IIIIl we select several model points and estimate the 
expected number of the stopped particles for J L = 300 fb _1 and 3000 fb _1 . In Section HVT 
we discuss the measurement of the two body decay f — > rX. The three body decay of the 
CNLSP is studied in Section |VJ Section |VTJ contains the discussion and comments. 
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diameter 


weight of the detector 


length 


ATLAS 
CMS 


22m 
15m 


7Kt 
12.5Kt 


44m 
21m 



II. ASSUMPTIONS ON THE STOPPERS 

In this section, we discuss the possibility to install massive stoppers next to the LHC 
detectors, CMS and ATLAS. It turns out that the CMS cavern may allow for an easier 
installation and more room for a massive stopper, compared to the ATLAS cavern. The 
parameters of the two detectors are listed in Table. [Ill The diameter of the CMS detector is 
smaller, therefore the massive stopper can be placed closer to the interaction point at CMS. 
The weight of the CMS detector is about twice as large as that of ATLAS. Because of the 
large weight, the cavern of CMS is designed so that it can sustain a massive object safely, 
which includes a reinforced floor to spread the detector pressure more equally. A potential 
massive stopper with a weight of around a few kton can be placed on both sides of the CMS 
detector, but it will need a reorganization of the scaffolding and gallery paths on the cavern 
walls, to make room for such an additional detector. 

The assembly and construction of the two detectors is also very different. Most of the 
CMS detector components are assembled on the ground, and about 15 large detector units 
will be lowered in the cavern for final assembly works. Hence the installation is less integrated 
with the cavern, leaving relatively more freedom and thus changes needed in the cavern to 
install massive stoppers are somewhat simpler Q. On the other hand ATLAS detector is 
assembled mostly in its cavern. The huge magnets toroids and outside muon system fill up 
the cavern. The cryogenic system in the ATLAS cavern is also taking space outside the 
detector 

We assume two stoppers with the size 3.5m x 15m x 15m and the average density 
Pstop = 5g/cm 3 , hence, the total weight of the detector is 8 kton. This is maximum possible 
rectangular parallelepiped volume that can be placed in the cavern, with the given space 
to the cavern wall and with its long edge being the same as that of the CMS barrel part 
(see Fig. [p. The stoppers are thus placed 8.5 m away from the interaction point. We also 
assume that stopping power of the CMS detector is equivalent to (2500/ sin^) g/cm 2 iron, 
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FIG. 1: Left: a schematic figure of the CMS detector and two stoppers. The numbers are in units 
of meters, and (0, 0, 0) is the collision point. Right: two stopper-detectors and a circle about the 
size of CMS detector are superimposed on the cross section of CMS cavern UXC 55, drawing taken 
from Ref. Q]. 

where 9 is angle between the CNLSP direction and the beam direction. The number comes 
from the average density of CMS detector, 3.37g/cm 3 , which leads to the weight per cm 2 
for the radial direction of 2500g/cm 2 . 

As discussed in the previous paper [J] , the stopper can be a hadronic and electromagnetic 
calorimeter simultaneously, if the detector consists of layers of dense stopper and tracking 
devices. The measurement of the energy of the decay product of the CNLSP is the key 
ingredient to explore the CNLSP interactions to the X particle. In this paper we assume 
that the CNLSP is the scalar tau lepton f, which decays mostly as f — > tX where X = G 
or a. 

The r decays into Wii> T , or into n ± and 7r°'s. We do not consider the decays into /i, 
because the muon energy cannot be measured unless the stopper contains a magnetic field. 
The energy of the leptons are much softer than the parent r energy anyway, so that they 
are less useful for the study of the decay kinematics. 

A large volume detector is advantageous to measure the energy of the r decay products, 
because the detector must contain most of the energy of the showers from the r decay prod- 
ucts. To fully absorb the hadronic cascade one needs sufficient thickness of the calorimeter. 
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The required thickness for an iron calorimeter (density is 7.87g/cm 3 ) is listed in Ref. 101. 
and is about 170 (120) cm for 100 GeV single hadron energy for 99% (95%) containment 
respectively, equivalent to 1337.9g/cm 2 (944.4g/cm 2 ). To measure the energy with suffi- 
cient accuracy by a stopper with p = 5g/cm 3 , the distance between the decay position 
and the end of the detector along the shower path must be at least 190cm. A simplified 
Monte Carlo simulation shows that 63% of the f decays satisfy this condition. High energy 
photons and electrons initiate electromagnetic cascades, which are much easier to contain 
in a detector. The energy deposition along the axis of the cascade is well described by 
a gamma distribution, where the maximum occurs at ~ 8Xo, where Xq is the radiation 
length. The energy deposition terminates around 15X . (See Fig. 27.17 of 10J.) For the 
case of iron (X = 13.84g/cm 2 ), 200g/cm 2 thickness of the material is needed to stop the 
electromagnetic cascade. The difference in the absorption length is used to discriminate 
isolated photons from hadron, which would be useful to study f — > r^X decay. 

In a previous paper we discussed the possibility to re-use of the existing 1 kton 
detector, such as SOUDAN II [l2| as the CNLSP stopper. The SOUDAN II consists of the 
layers of 0(m) long drift tubes and thin iron plates. The physics goal of this detector is 
the search for proton decays. To be sensitive to the low energies involved in these decay 
processes, the average density of SOUDAN II is low, less than 2g/cm 3 . The size of the 
detector is probably enough to stop a certain amount of CNLSPs and measure the decay 
rate, and therefore it may well be appropriate for the a first stage of the CNLSP study. On 
the other hand, it is certainly not enough for a detailed energy measurement of the f decay 
product. Most of the hadronic decay cascades are not fully contained for a detector with 
the geometry allowed by the space in the CMS cavern. A high density detector consisting 
of the layers of drift tubes/scintillators or RPCs to measure the charged particles between 
iron plates thicker than SOUDAN II will be more optimized for the CNLSP study. 

In this paper we assume a conservative energy resolution for hadronic showers, which is 
around 150%/ \jEj GeV. The value is not unrealistic for a simple massive and affordable 
detector, if a shower is sufficiently contained in the stopper. An additional complexity 
may occur for showers which develop parallel to the layers of the tracking devices. If the 
particles pass mostly through tracking devices, they feel a much lower average density, while 
the particles going mostly through the iron plates do not give detectable signature efficiently. 
We assume that the measured energy will be corrected depending on the shower directions. 
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A detector uniform to all directions would be better to measure the energy of decay products 
of the stopped CNLSP. The calorimetery technologies studied for the ILC may satisfy such 
conditions; see [3] for the CNLSP study at ILC. 



III. SUPERSYMMETRIC MODELS WITH CHARGED NEXT LIGHTEST SUSY 
PARTICLES, AND EXPECTED NUMBER OF STOPPED CNLSP 

In this section, we briefly describe supersymmetric models with a long-lived charged next 
lightest SUSY particle. We also select some model points, and estimate at each model point 
the number of CNLSP that can be stopped at the stopper-detector. 

In the minimal supergravity (mSUGRA) models, the scalar masses and gaugino masses, 
trilinear couplings are universal at the GUT scale Mqut which are denoted by mo, M1/2, 
and Aq respectively. The resulting mass spectrum is obtained by solving renormalization 
group equations (RGE's). When the RGE is integrated up to 0(m z ), slepton soft masses 
and gaugino masses are approximately given by the following convenient formulas 

m\ L ~ ml + Q.3M* /2 , 
ml ~ m~ ~ ml + 5.4Mf 



mf ~ m 2 + 0.5Mjy 2 , 
™h ~ ml + 0.15Mf /2 , 



e R ~ " fr ' u - L ^" J l/2) 

Mi = Mm m 

9 2 ' \ / 

9i 9x 

where gx is the gauge couplings at the unification scale. 1 In this model, the mass of gravitino 
m^ is order of m or M1/2, because both of them are proportional to the F /M p \ with 0(1) 
coefficient, where F is the fundamental SUSY breaking scale. Depending on the 0(1) 
coefficients, the gravitino can be the lightest superparticle in general. Also in the gaugino 
mediation models [14], the gravitino can be the LSP in a large domain of the parameter 
space [l^]. In this case the scalar masses are very small at the boundary, i.e., mo — 0, and 
the stau naturally becomes the NLSP. 

In the gauge mediation (GM) model [3] the supersymmetry is broken at lower energy 
scale, and the SUSY breaking in the hidden sector is mediated to the MSSM sector by gauge 



1 The pole masses of strongly interacting SUSY particles receive a large corrections of 0(30%) if the mass 
scale is O(lTeV). 
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interactions. The simplest GM model is described by 6 parameters: A = F/M mes , M mes , N 5 , 
tan/3, sgn/j = ±1 and gravitino mass. Here iV 5 is an effective number of vector-like heavy 
quarks and leptons in SU(5) representations, $ and $. (N 5 = 1 for 5 + 5 quarks and leptons, 
and N 5 = 3 for 10 + 10.) A messenger field Y couples to the vector-like fields as W = XY 
and develops a vacuum expectation value X(Y) = M mcs + OOF. The gravitino mass is given 
by rriQ = F /(y/3M p i) where F (> F) is the total SUSY breaking scale of the theory. We 
take rriQ as a free parameter in this model. 

The masses of MSSM particles at M mes are obtained by relatively simple formula. Gaugino 
masses satisfy the GUT relation, and 



m-g 



(2) 



Squark and slepton masses are given by 



«4 



m ln 



m 2 i 

an 



m e R = 



3 V4tJ 



2 + 3 



2 V4tt) 
8 /a s y 8 fai\ 2 ' 
3 \4n) + 15 \4tt) 

3 Utt/ 15 Utt/ 

2 Utt/ 10 Utt/ 
^^ 2 liV 5 A 2 . 



+ 



J_ /«iV 
30 Utt/ 

iV 5 A 2 , 



AT.A 2 , 



iV 5 A 2 , 



^,A 2 , 



(3) 



5 V47T/ 

All of the above models predict a large mass hierarchy between strongly interacting 
superpartners and weakly interacting superpartners. Heavy gluino and squarks are copiously 
produced at the LHC, and they decay into the light weakly interacting SUSY particles. The 
lightest SUSY particle in the MSSM sector is either the lighter stau fi or the lightest 
neutralino Xi- 

The mass of the fi is the smaller eigenvalue of the mass matrix, 



+ m 2 - 1(277^ - m|) cos 2/3 
—m T (A T + fi tan j3) m 



—m T (A T + /j.tan/3) 
+ m 2 + [mly — m|) cos 2(3 



(4) 



Because of the off-diagonal elements of f mass matrix, f-y could be significantly lighter than 
the other sleptons. If the f x is the NLSP, the stopper-detector is useful to stop it and to 
study its decay. We therefore consider the phenomenology when m g ,mg ^> m^o > in 
this paper. In the following, we omit the subscript of fi and denote the NLSP stau as f . 



In most part of this paper, we will discuss the CNLSP physics as model independent as 
possible. However, to give some numerical reference, we choose a few model points. For 
mSUGRA models, we take the points proposed in [4J. Theparameters are listed in Table. [U 
and mass spectrum is given by the ISA JET ver. 7.69. In [4j, the gravitino mass is taken to 
be tuq = m in those model points. 
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m^o [GeV] 


175 


417 


416 


m h [GeV] 


154.2 


343.5 


324.3 


cr(SUSY)[pb] 


3.03 


2.27 x 10" 


" 2 2.34 x 10~ 2 


stopped in the stopper-detector per 10 5 events 


255 


250 


254 


stopped for 300 fbr 1 


4636 


34 


36 



TABLE I: Some model points in mSUGRA model from [4|. The mass spectrum and production 
cross section relevant to our study are shown. 



A[TeV] 


40 50 60 70 80 


m~ g [TeV} 


0.93 1.13 1.34 1.54 1.74 


m^o [GeV] 


161.7 205.3 248.7 292.1 335.4 


m fl [GeV] 


120.5 150.1 179.9 209.8 239.9 


a(SUSY)[pb] 


5.24 1.68 0.64 0.28 0.13 


stopped in the stopper-detector per 10 5 events 


282 274 274 294 302 


stopped for 3001b"" 1 


8830 2762 1052 494 236 



TABLE II: Some model points in gauge mediation model. The production cross section and mass 
spectrum relevant in this study are also shown. 

For the gauge mediation model, we take the model points similar to that for the 
study in [l|. Namely, we fix N 5 = 3, tan/? = 15 and A/M mes = 0.5, where A = 
40, 50, 60, 70, 80 TeV. The mass spectrum and production cross sections are summarized 



in Table. [TTJ The gravitino mass for these model can be very small, whose minimum value is 
given by 0(A 2 )/M p \. If the messenger sector couples a fraction of total SUSY breaking, i.e. 
F < Fq, the gravitino mass, which is proportional total SUSY breaking F , can be large. 

The pattern of the mass spectra are similar for both mSUGRA and GM sample points. 
The gaugino masses obey the GUT relation, and f mass is lighter than B mass because of 
the CNLSP assumption. For mSUGRA, this implies mo < M1/2, so that q mass is about as 
large clS Q ? clS it is so in GM models. The relation rriq ~ m- g ^> m? is realized in our model 
points. 

To estimate the number of stopped CNLSP, the production cross section of the SUSY 
particles and the velocity distribution of f at LHC must be evaluated. They depend on the g 
and q masses, and the mass difference rrig^ — rut- We estimate the production cross section 
and the decay distribution by using HERWIG[17], where the mass spectrum and branching 
ratios are interfaced from ISAJET[18] to HERWIG by using ISAWIG[19]. We generated 10 5 
SUSY events for each model point. 

The flying range R of the charged stable massive particle may be calculated by integrating 
the energy loss equation of heavy ionizing particle (Bethe-Bloch equation). The result is a 
function of (3 = p/E, with a linear dependence on the mass of the particle M. In this paper, 
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See also the 



the stopping power of the stopper-detector is calculated using the data in 
detailed discussion in Ref. 3]]. 

The maximum length of a particle track through the stopper-detector /(max), which 
depends on CNLSP direction, is calculated assuming that the track of the CNLSP goes 
straight from the production point. We regard the particle is stopped in the stopper-detector 
if the flying range R satisfies: 

2500 g/cm 2 n 2500 g/cm 2 

r^- < R < ^- + / max x p stop , 5 

sin sin 

where p st0 p = 5g/cm 3 is the density of the stopper-detector. The number of stopped CNLSP 
in a stopper per 10 5 events, and the number of stopped CNLSP for £ = 300fb _1 in the two 
stoppers, are listed in Table HI and HT1 In Fig. [2] we also show the simulated stopped positions 
in a stopper-detector. The position distribution is uniform in the detector. 

In Fig. [31 we show the /?7 distribution of the CNLSPs for a few model points. We find 
that CNLSP tends to be less relativistic when the SUSY scale is large, because momentum of 
parent squarks and gluino is reduced. In Fig. [3], the peak position of the is at ~ 1.5(1.3) 

10 
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FIG. 2: The simulated positions of the CNLSPs in a stopper-detector for the A = 40 TeV GM point. 
Here the z axis is the beam direction, and the y axis is the vertical direction. The origin (0, 0, 0) is 
the interaction point and we assume the stopper-detector is at 8.5m< x <12m, — 7.5m< y <7.5m, 
and — 7.5m< z <7.5m (cf. Fig. [T]). The big square is the projection on the y-z plane, the top 
rectangle is the projection on the x-z plane, and the right rectangle is the projection on the x-y 
plane. 

for A = 50(80) TeV. Because of that, the number of CNLSP in the smaller /3"f region is 
increased as the SUSY scale is increased. On the other hand, as mcNSLP increases, CNLSPs 
with smaller /?7 are stopped in the detector while the number of events in the smaller (3 r ) 
bins are kinetically suppressed. For instance, for A = 50 and 80 TeV, CNLSPs in the 
bins between 0.5 < (3-y < 0.6 and 0.45 < fiy < 0.55 are stopped in the stopper-detector, 
respectively. Altogether, the number of stopped f for A = 40 to 80 TeV for 10 5 SUSY events 
is roughly constant as we can see in Table. HT1 

The production cross section reduces when the gluino mass is increased because the parton 
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FIG. 3: The f3j distribution of f for A = 80 TeV and A = 50 TeV in GM models and for point C 
in mSUGRA models for 10 5 SUSY events. Light (dark, gray) shaded histograms are the number of 
stopped events in the stopper scaled by factor of 10, for point £ and for GM models with A = 80, 50 
TeV respectively. 

distribution functions of gluon and quarks are quite small for x ^> 0.1. For 300fb _1 , the 
number of events stopped at the assumed two 4 kton stoppers ranges from 8000 events to 30 
events in the table. We will see in the next section that accumulation of (9(100) CNLSPs are 
enough to measure the lifetime with (9(10)% accuracy. We will also estimate the resolution 
of r lepton energy E T arising from the decay f — ► rX through the end point of the tau 
jet energy. Statistically 0(1000) stopped CNLSPs are enough to measure the end point 
with a few GeV error. In the upgrade of LHC (SLHC), integrated luminosity of 3000fb _1 is 
proposed, therefore the number of stopped events ranges from (9(300) to (9(80000) for the 
model points presented in Table [Til 

IV. STUDY OF THE f TWO BODY DECAY IN STOPPER DETECTOR 

In this section we study the two body decay of the CNLSP in the stopper-detector. Both 
in mSUGRA and GM models, the stau can be the CNLSP and decays into the gravitino G 
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FIG. 4: The lifetime of the CNLSP f in the case of gravitino LSP, as a function of the gravitino 
mass, for the stau mass 100, 200, and 400 GeV (from top to bottom). 

and the T-lepton. The CNLSP decay width into a gravitino and a lepton is given by 
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T 



rnQ + m 2 



48ttto|M 2 



m- 



(m? — m?£, — m 2 ) 2 
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m 2 £ + m 2 
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m G J 



[m- 



ml 



2\2 



3/2 



(6) 



We show the dependence of the stau lifetime on the gravitino mass in Fig. HI 

In general, the two-body decay f — > Xt can be triggered by a single tau jets initiating 
from the position where f is stopped. The tau energy is monochromatic and expressed as 



m\ + m 2 . 



rn 



x 



2rrii 



(7) 



Here X is the invisible particle in the f decay, in this case X = G. If one can measure both 
of the lifetime and the mass of the stau, the gravitino mass can be determined assuming that 
Eq.fJH]) is correct. Then the total SUSY breaking scale F = x/Sm^Mpi is also determined, 
which is very important to understand the hidden sector physics. 

The LHC main detectors can determine the mass of f through the stau velocity mea- 
surement (3f, e.g., in the muon system of the CMS detector. However, the measurement of 
the lifetime may not be easy at the main detectors if the lifetime is too much longer than 
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the detector size. The cross section is typically 0(1) pb or less, therefore we have at most 
O(10 5 ) event at hand for C = lOOfb -1 . Thus, only 10 events or less decay inside the detector 
if cr > 100 km. Some of the CNLSP are stopped in the main detector, but measuring the 
decay precisely in the main detector would be challenging during the beam time. On the 
other hand a stopper-detector [2j can measure the position and the time where a CNLSP 
is stopped, and its decays without dead-time. The lifetime measurement will be discussed 
in Sec JIV Al See [4] also on the idea to measure the lifetime by triggering muons from the 
decays of the CNLSP stopped in the surrounding rock. 

To predict the CNLSP lifetime, one has to determine the gravitino mass independently. 
This is possible through the extraction of E T from the energy distribution of the tau jet 
from the CNLSP decay, because E T is a function of m? and mx, as can be seen in Eq. (j7j). 
mx is expressed function of E T as follows; 



mx = y rri~ + m 2 — 2m^E T . (8) 

Fig. [5] shows the dependence of the reconstructed LSP mass Eq.(jHJ) on the tau energy E T 
for several values of m f . As can be seen from the figure, it is crucial to measure the E T 
as precise as possible, especially for small mx, in order to determine the mass mx- Hence, 
the stopper-detector should offer a reliable measurement of E T . We will discuss the E T 
measurement and the mx reconstruction in Sec. IIVBI and Sec. IIV C\ respectively. 

One can study the supersymmetric version of gravity interaction by studying the consis- 
tency between the observation and prediction of the decay rate. If the gravitino is the LSP, 
the decay width is given by Eq.(jSJ). Now if one can independently determine the gravitino 
mass by means of kinematics, as described above, Eq.© can be used in the other way round, 



which leads to the measurement of the 'supergravity Planck scale' j^] 



ml ( _ m 2 d + m 2 T 



M pl (supergravity) = — — ^— 1 ^—^ 



4 r a 2 2 1 3 / 2 

Am-m T 

1 ' 



(mf — TOg, — m 2 ) 2 



(9) 



Comparison of the obtained M p \ (supergravity) with the Planck scale of the Einstein gravity 
M pl (gravity) = (SttG)- 1 / 2 = 2.43534(18) x 10 18 GeV Q would be a crucial test of the 
supergravity. Prospects of the "Planck scale" measurement will be discussed in Sec lIVDl 

It should be noted that the undetectable particle X may not be the gravitino G. Any 
particle which couples weakly to f can be particle X. If the decay width into Xr is larger 
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FIG. 5: Sensitivity of the reconstructed LSP mass mx [see Eq.([8])] on the tau energy E T , for 
Tfif = 100, 150, 200, 250, and 300 GeV (solid lines). Dashed lines show contours of the lifetime in 
the case of the gravitino LSP (X = G). 

than the decay width into Gt, i.e. T(f — > rX) > r(f — > tG), 2 the CNLSP lifetime may be 
different from the supergravity prediction obtained from measured m f and tau jet energy 
distribution. Inconsistency between the measured and predicted lifetime immediately means 
a discovery of a new sector that may not be accessible otherwise. 

One of well motivated examples of such a non-SUGRA decay is f — > ar where a is the 

n 

axino, superpartner of the axion. The CNLSP r decay into axino is studied in Ref. [6j for 
hadronic, or KSVZ axion models 21]. In this paper we adopt the set-up in Ref. jf| for the 



axino interaction, which we briefly describe here. 



In KSVZ axion model, the Peccei-Quinn (PQ) mechanism [22] is realized in an action 
with new heavy quarks. When these heavy quarks are integrated out, anomalous terms 
involving the axion and gauge bosons are generated at low energy effective action. When 
the axion interaction is supersymmetrized, its fermionic superpartner, the axino a, must be 
introduced. The axino mass m„ can range between the eV and GeV scale depending on 



the model and SUSY breaking scheme 23j, and we regard it as a free parameter in this 



paper. The coupling of the axino to the bino and the photon/Z-boson at the scale below 



2 This includes the trivial case where the decay into the gravitino is kinetically forbidden, T(f — > tG) = 0, 
i.e., m,Q > mf > mx- 
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FIG. 6: Feynman diagrams for f — ► ar. 
the Peccei-Quinn scale f a is given by the Lagrangian 

L ~a = i i* a ^2 m ^idl n^AB (cos W F^ - sm9 w Z^). (10) 

167T COS Ofyfa 

The action does not contain direct fra coupling and also strongly suppressed by the PQ 
scale 10 9 GeV < f a < 10 11 GeV. The two body decay f — > ar is induced by the one loop 
diagram shown in Figj6l The loop integral has a logarithmic divergence. This is because the 
effective vertex ([TO]) is applicable only if the momentum is smaller than the heavy (s) quark 
masses, whereas the loop momentum exceeds that scale. In the full theory, one calculates 



two-loop diagrams with the heavy (s)quarks, which leads to a finite result 24(. Here we 



regulate the logarithmic divergence with the cut-off ~ f a p, |25(], and the effective arfn 
coupling is parameterized as 

L ~ -£C aYY — 2 tw~T lo § —TRTP L a + he, 11 

StHcos 4 ^ f a m 

Here, m ~ m~ g ~ (9(100 GeV) and we take \og(f a /m) = 20.7 hereafter. The parameter £ 
is an order one parameter to represent the uncertainty coming from the cut-off procedure 
mentioned above, 
width is given as jg] 



n this paper, we regard this as a free parameter. The two body decay 

El 



1 ' 512^ cos 8 ^ fl m% ^ g \m 

_ « 25 secr^ y (l - (^) (^) 2 (^) 2 



16 




FIG. 7: An example of maximum-likelihood fit of the temporal distribution N(t) for N tota \ = 100. 
In the left plot, the solid line is the generated events, and dashed line is the best fit. 

A. Lifetime measurement 

In this subsection, we estimate the statistical error in the CNLSP lifetime measured by 
the stopper-detector. The analysis is model-independent and applicable to any long-lived 
charged particle stopped in the stopper. 

For each CNLSP stopped in the stopper-detector, the stopping time t stop and the de- 
caying time tdccay will be recorded. The lifetime is then measured by fitting the temporal 
distribution of the decaying events N(t), where t = tdecay — t stop . Here, we use a maximum- 
likelihood fitting and adopt the following procedure: First we define a time t e such that 
N(t < t e ) = (1 - e" 1 )^! ~ 0.632AT total and N(t > t e ) = e^A^i - 0.368AT total , where 
-^Vtotai is total number of stopped event. 3 (For large N tota \, this t e is already a good estimator 
of the lifetime.) We then calculate the InL distribution as a function of lifetime r: 

In L(t) = fp(jk-Mr)) (13) 

i=bins 

where /p(n; v) = v n e~ v jn\ is the Poisson distribution, rii is the number of events in the i-th 
bin, and Viij) is the predicted average number of events in that bin: 

Viir) = TVtotai (e" (i - 1)A * /T - e- iA ^) ■ (14) 
Here, we take 1 bin = At = t e /5. 

3 More precisely, t e is defined by t e = (tj + tj + \)/2, where j < (1 — e _1 )iVtotai < j + 1 and tj is the decay 
time of the j-th event: t\ < ti < ■ ■ ■ < tj < ■ ■ ■ < £jv total - 
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An example of the time distribution and the best fit curve for MC distribution are shown 
in Fig.[7]for A^otai = 100. The no confidence interval can be estimated by the range 2 A In L = 
2(lnL max — InL) < n 2 . Fig. [7] shows that the error of the lifetime is Ar/r = (10 — 15)% for 
-Wtotai — 100. We have done the same analysis for iV total = 1000 and found Ar/r = (3 — 4)%. 

So far, we assumed that most of the stopped CNLSPs decay within the experimental time 
scale. We now consider the case where the lifetime is longer. Suppose that one observes 
IVtotal = 1000 stopped CNLSPs and only 10 events decaying within 1 year. In such a 
case, the lifetime is estimated from the number of the decaying events. For instance, 95% 



interval of the mean v of Poisson variable for n = 10 is v = [5.4,17.0] [10] . Using v = 
Atotai(l — e" 1 y ear / T ), a 95% interval 58 < r < 184 years is obtained. For much longer 
lifetime or smaller iV tota i, only a lower bound of the lifetime is obtained. 

We have assumed that the background is negligible. The background from cosmic neu- 
trino and hard neutrino produced from the main detector interaction point is small (cf. j^j]), 
however careful study on the accidental background is necessary when statistics is low. This 
is beyond of the scope of this paper. 

B. Measurement of the r energy from distribution of the r jet energies 

In this subsection, we estimate the uncertainty of the tau energy determination. Then, 
in Sec JIV Cl we discuss the kinematical reconstruction of the LSP mass. Schematically, the 
procedure is as follows: 

.Ejet distribution — > E T — > m X - (15) 

When f decays into r and invisible particle X, the tau energy E T is monochromatic (see 
Eq.flTj)). E T can be obtained by fitting the r jet energy distribution : 

jyy decay products 

(£ jet ; E T ) where E jet = V E t . (16) 



dE ]ct K Jct ' T! Jet ~ 

Among the decay products of the r lepton we omit the neutrinos and muons. In order to 
see the prospects of E T measurement with a finite number of events, we generate the events 



from r decay by using the TAUOLA 



261 ]. and we perform a maximum-likelihood fitting of 



low statistics ("experimental") event sets by high statistics ("theoretical") distributions. 
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FIG. 8: Tau jet energy distributions for left-polarized (left) and right -polarized (right) tau gener- 
ated by TAUOLA. The primary tau energy is E T = 50, 70, 90 GeV from left to right. The dashed 
lines show the spectrum without energy resolution effect, and the solid lines show the spectrum 
with detector energy resolution /S.E- ]Ct / E- ]Ct = 150%/ £j ct /GeV. 

In practice, the observed jet energy distribution depends not only on the E T , but also on 
the detector resolution AE^ et and the tau polarization P T : 

-^-(£ jet ; E T ,P T ,AE jet ). (17) 

To obtain theoretical predictions for dN j ' dE- iCt (E- iCt \ E T , P T , AEj et ), we have generated tau 
decay events with high statistics run of TAUOLA for the parameter space E T = 30 - 
125 GeV and -1 < P T < 1, while fixing AE ]ct /E ]ct = 150% /^/E~jGeV 4 For each single 
t decay generated by TAUOLA, the jet energy E- ]et is calculated, and then smeared by a 
Gaussian fluctuation with a variance a 2 = (AEj et ) 2 . Fig [H] shows examples of the jet energy 
distribution. 

In order to see a realistic situation, we generate a small number of tau lepton decays 
for a fixed parameter set of (E T , P T , A£j et ), and then fit the result by the "theoretical" 
distribution obtained above. Here and hereafter, we assume that the energy resolution 
will be known in advance. Examples are shown in Fig. [9] and Fig. [10] for total number of 
decaying tau leptons N T = 1000 and A^ T = 100, respectively. The energy of r jets which 
are not contained in the detector may not be measured precisely. The N T corresponds to 

4 We have generated 10 6 events for each of the parameter sets E T = 30, 31, 32,... 109, 125 GeV and P T = ±1 
(i.e., 96 x 2 parameter sets x 10 6 events), and interpolated the distribution between those parameter points. 
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FIG. 9: An example of maximum-likelihood fit of low statistics events from tau decays, (a) 
Energy distribution of tau jet events generated from 1000 decaying tau leptons for E T = 70 GeV, 
AEjet/Ejet = 150%/ \/ Ej ct I GeV and P T = +0.8 (solid histogram), and the best fit distribution 
(dotted histogram). Only the bins between the vertical lines are used for the fit (see text), (b) 
2AlnL = 2(lnL max — InL) projected onto the E T axis, (c) Contour plots of 2AlnL = 1, 4, 9 
projected onto the (P T , E T ) plane. 

the number of well contained events. Some comments are in order here, (i) We take 1 
bin = 10 GeV. (ii) For the maximum-likelihood fitting we take only the bins with number 
of events > 1, and the bins above the peak energy, which in the examples of Fig. E(a) and 
HPT a) correspond to the bins between the two vertical lines, (iii) We then calculate the In L 
distribution in the parameter space of (E T , P T ) as follows 

In L(E T , P T ) = In /P (V W ;Nt&(E T ,P T j) (18) 

j=bins 

where /p(n; v) = u n e~ u /n\ is the Poisson distribution, Nj™ is the number of events in the 
z-th bin for the low statistics run, and N^ lgh (E T , P T ) is the predicted number of events in 
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FIG. 10: An example of maximum-likelihood fit of low statistics events from tau decays. The same 
as Fig. [9] but with 100 decaying tau leptons. 

the i— th bin as a function of the parameter set (E T , P T ), normalized by the total number of 
decaying tau leptons. Here, we assume that the total number of stopped f is known, i.e., 
the total number of events is not taken as a free parameter for the fit. 

The no confidence interval can be estimated by the range 2 A InL = 2(lnL max — In L) < n 2 
projected onto the E T axis [see Figs. MJo) and ITUT b)]. One can see that the primary tau 
energy can be determined within an error of a few GeV. From Fig. M^c) one can also see 
that the polarization is hardly determined by the energy distribution analysis even with 
N T = 1000. This is because the sensitivity to the polarization becomes very weak once the 
finite detector resolution is taken into account (cf. Fig. [8]). 

To estimate the statistical error in E T measurement we have generated the event sets with 
same statistics N T = 100 (1000) and repeat the fit on E T and P T to obtain the best fit value 
of the E T , E* est . In figure HU we show the distribution of E* cst for E T = 50,70,90 GeV, 
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FIG. 11: Distribution of E^ est for E T = 50, 70 and 90 GeV (from left to the right), for N T = 1000 
(left) and N T = 100 (right). 



P T = 1, and A£ jet /£ jet = 150%/^/E/GeV. From the variances of these distributions, we 
estimate the la uncertainty of E T as 



5E T /E T = 61%/V^r/GeV for N = 100 , 



6E T /E T = 15%/V^r/GeV for = 1000 . 



(19) 
(20) 



C. Determination of the LSP mass 



Now we can estimate the uncertainty of the LSP mass. In Fig. [12] we plot the range of 
reconstructed LSP mass 



m x 



rrif 2 + ml — 1m^E T 



(21) 



where 



E r -5E T (E T ) < E T < E T + SE T (E T ) 



(22) 



and E T = (m^ + m^. —m 2 x )/2mf. Note that the reconstructed LSP mass depends on not only 
E T but also the measured CNLSP mass m~f. In Fig. [121 we show the range of LSP mass for 
mf = rrir (solid lines) and 0.99mf < rrif < 1.0 lm? (dashed lines), the latter corresponding 
to 1% uncertainty of the CNLSP mass. The stau mass determination from time of flight 
was discussed in |l|, |4|, [27J. For CMS detector the mass resolution is estimated as 10-20% in 
each event and less than 1% for ~ 1000 events. 
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FIG. 12: Uncertainty of the reconstructed LSP mass, corresponding to the estimated la uncertainty 
of E T , for N T = 1000 and 100, and mf = 150,200 and 300 GeV. Solid lines represent the case 
without an error of mf, and dashed lines include 1% uncertainty of mf. Vertical dotted lines 
represent the stau lifetime of 1 month, 1 year, and 10 years in the case of gravitino LSP {X = G). 

As can be seen in Fig. [121 the kinematical reconstruction of the LSP mass is possible 
if mx is sufficiently large, mx ~ 0.15mf for Nf = 1000 and mx ~ 0.25mf for N? = 100. 
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Otherwise one can get only an upper bound on the mass mx- 



D. measurement of the "Planck scale" 

Finally, if the LSP is the gravitino (X = G), the uncertainty of the reconstructed gravitino 
mass translates into an uncertainty of the supergravity Planck scale, which is obtained 
by substituting for tUq in Eq.([n]). To take into account the error of m? , we also substitute 
mf for f)if. Eq.([9]) then becomes 

2 1 (m f E T - ml) (e t 2 - m 2 T ) 1 

M - = ^ ~,+ ,,~r • (23) 

< 37ri t mf + m l T — 2m?rj T 

which is shown in Fig. [13] for m~f = mf (solid lines) and 0.99mf < mf < l.Olmf (dashed 
lines). Here we do not include the uncertainty of the lifetime measurement, which simply 
affects the measured Planck scale as M p i oc (l/Tf) 1 / 2 . Neglecting the r-lepton mass, Eq. (123]) 
is simplified as 

- — -4 

gi 2 * , * ^ ■ (24) 
3nr f ffe _ 2E T 

As discussed in the previous section, the kinematical reconstruction of the gravitino mass 
is possible only if mg is sufficiently large, > 0.15mf for N? = 1000 and m^ > 0.25m ? for 
Nf = 100. For smaller values of the gravitino mass, one can get only a lower bound on the 
Planck scale M pl . 

One can see whether the determined 'Planck scale' is inconsistent with the Planck scale of 
the Einstein gravity, M p \ = 2.4 x 10 18 GeV. In other words, one can test the assumption of the 
decay f — > tG by comparing the observed lifetime with the predicted lifetime. For example, 
if the NLSP dominantly decays into axino, the lifetime would be far shorter compared with 
gravitino assumption for a fixed mass mx, leading to a smaller value of fitted M p \. Suppose 
for instance one measures m f = 150 GeV and mx = 30 GeV. When X — a, the lifetime 
becomes O(10 sec) for f a ~ 10 11 GeV and m^ ~ m f [cf. Eq. (fT2l) ]. If one uses Eq. (1231 . 
a "Planck scale" M p i = O(10 15 GeV) would be obtained, thereby falsifying the gravitino 
assumption. 
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FIG. 13: Uncertainty of the Planck scale measurement, corresponding to the estimated la uncer- 
tainty of E T , for N T = 1000 and 100, and m f = 150, 200, and 300 GeV. Solid lines represent the 
case without an error of mf measurement, and dashed lines include 1% uncertainty of m,f . Vertical 
dotted lines represent the stau lifetime of 1 month, 1 year, and 10 years. 



25 



E. model points and cosmological constraints 



We now discuss mSUGRA model points discussed in section UTTj together with cosmolog- 
ical constraints. In the early universe, the f CNLSP has been in thermal equilibrium until 
its decoupling, T d ~ 0.04m f . If such a particle has a very long lifetime, as discussed in this 
paper, its decay during or after the big-bang nucleosynthesis (BBN), Tbbn ~ 1 MeV, may 



spoil the successful prediction of BBN 
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29|. In the models with f NLSP and G LSP, 



this leads to severe constraints on the parameter space of fmf, mg), in particular to upper 
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bounds on the gravitino mass for a given stau mass 

Furthermore, it has been recently pointed out that a heavy charged particle can form 
a bound state with light elements during the BBN, which can lead to new effects and/or 
severer constraints [32j. In this paper we do not discuss these effects because it is difficult 
to evaluate the net effect quantitatively and it still awaits detailed analysis. 

We should also mention that those BBN constraints may disappear if there is entropy 
production between the stau decoupling (Td ~ mf/20)_and the BBN (Tbbn ~ 1 MeV), 
because the stau abundance is diluted before its decay [331 ] . 

Keeping in mind the possibilities of severer bounds and also a possible loophole, let us 
discuss the cases of nif = 150 GeV and nif = 340 GeV, corresponding to the mSUGRA 



model points e and £. According to the latest analyses [3l| including the effects of the 
hadronic decay [29], the bounds on the gravitino are tuq < (20-80) GeV for ~ 150 GeV 
and Trig ~ (40-200) GeV for nif ~ 340 GeV. 5 The ranges of upper bounds correspond to 
the uncertainties of various bounds from primordial light elements. 

For rrif = 150 GeV (model point e), as can be seen in Tabled! one could collect more than 
1000 staus for 300 fb^ 1 . The bound txiq < (20-80) GeV then suggests that the measurement 
of the gravitino mass and the "Planck scale" may become possible if one assumes conservative 
BBN bounds and if the gravitino mass is sufficiently large (cf. Figs. [12] and [T3|) . For 
rrtf = 340 GeV (model point (), the measurement would become easier if one could collect 
the same number of staus (cf. Figs. fT2l and [TBI . However, from Tabled we find that expected 
number of stopped CNLSP is around 30 for 300 fb _1 . This is because gluino mass is above 2 
TeV for this point and production cross section is small. One needs SLHC (f £ = lOOOfb -1 ) 



5 These constraints were derived without using the bound on the 3 He. If one adopts it, the constraints 
become severer (cf. [29]). 
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to collect 0(100) events. 

In the case of axino LSP, the BBN bound is much weaker because the lifetime of the 
CNLSP stau becomes much shorter [see Eg. (1121) ]. Hence, the axino mass measurement is 
plausible for sufficiently large m„/m^. 



V. LIGHT AXINO VS GRAVITINO: THE RARE DECAY OF THE CNLSP 

A. Low energy effective action of the axino and gravitino involving photon 

When tuq ^ 0.2mf, it is difficult to determine the gravitino mass from r energy mea- 
surement at stopper-detector. Axino is a motivated candidate which couples weakly to 
the MSSM particle with comparable strength to the gravitino. In this section we therefore 
compare the decay f — > 0x7 with the decay f — > GV7. For simplicity, we will assume the 
NLSP is pure 'right-handed' stau, f = Tr throughout this section. Extension to the case 
with mixing with ft is straightforward, and the mixing angle dependence is expected to be 
small. 

The gravitino is a spin-3/2 particle. However, in the limit where tuq <C mf, the effective 
interaction to MSSM particles would be reduced to that of spin 1/2 particle, goldstino x- 
The effective action relevant to the tr decay is given as follows, 

2 

L = 7f \, $7"?* + H ^ + a m xb^YmcosOwF^ - ^9 W Z, V ). (25) 

V^m^Mpi 4v / 6m 3/2 Mpi 



The action is similar to that of axino given in Eqs. fllOp and ( flTl) except the coupling coeffi- 
cients. The relative weight of the two terms in Eq. (1251) are fixed by the supergravity, while 
for the axino the coefficient in Eq. (fTT]) is induced from Eq. fflQ|) by the radiative corrections. 
Note that the term proportional to j^BF^ is a non renormalizable coupling of the 

photon to gravitino or axino and induces significantly different 7, r distribution. 

The axino three body decay f — ► 7™ proceeds through the diagrams shown in Figfl^l 
where the hatched triangle express the effective vertex shown in Eq. ( TTTT) . On the other 
hand, the relevant diagrams for the three body decay into goldstino f — > 7r%, are given in 
FigJT5l The diagram corresponding to the top right of Fig. [TUdoes not exist for the goldstino 
case. The difference of the actions and the relevant diagrams will appear as the deviation 
of the decay distributions. In the Appendix, we list the three body decay differential width 
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FIG. 14: Leading Feynman diagrams for f — > rja. 



into gravitino/axino in the limit where the gravitino/axino mass can be neglected compared 
to ruf- The formula for the massive gravitino and axino are given in Ref. 6|. 



B. Numerical results 



The three body decay f — > t^X should be visible in the stopper-detector if it has an 
ability to measure charged tracks, and also segmented into small units. The position where 
f decays in the detector is the position where the charged track by the 7r + , ji and e from 
t is initiated. For hadronic tau decays, a tt^ is always in the decay products, sometimes 
with photons coming from tt° decays. Photons are converted into electron after passing ~ 
one radiation length Xq. Therefore for the iron based detector with p = 5(2)g/cm 3 , the 
photon shower starts 2.8(7)cm from the decay point. In summary, the three body decay 
of f —>■ tjX is identified as a charged track (which might be associated with collinear EM 
showers) + an isolated hard EM shower pointing back to the point where charged track is 
started. If the segmentation is not good enough, the efficiency to discriminate EM showers 
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FIG. 15: Feynman diagrams for f — ► tjG in the goldstino limit. 

from 7r ± would be reduced. 

As can be seen in the appendix, the three body decay amplitude can be written as a 
function of the angle between photon and r, 9, and E 1 . To be conservative we assume the 
energy resolution for the isolated photon shower is AE y /E 7 = 100% / \J E y / GeV and ignore 
the angular resolution of the photon momentum. (In the following E 1 denotes the photon 
energy after taking into account of this finite resolution effect.) Note that the shortest length 
of the detector is 3.5 m, which corresponds to 1750g/cm 2 for p = 5g/cm 3 . The EM showers 
are likely contained in the stopper because we only need 200g/cm 2 to fully absorb them. 

We need to require several cuts for the accepted events. 

• Experimentally, the angle between r and 7 must be large enough to avoid the overlap 
between r decay products and prompt 7. We only use the events where cos 9 < 
0.866. We do not lose sensitivity to the differece between the two decays f — > r^a/G 
by cutting these events. In the collinear region, the amplitude is dominated by the 
contributions from QED vertex, which is common for both of the decays. 

• The three body decay amplitude suffers soft and collinear singularity. Because we only 
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FIG. 16: The ratio R(X) = T(f — * Xjt)/T(t — » Xt) as a function of nig — nif. The solid, 
dashed, dot-dashed, short-dashed curves correspond to the case where nif = 100, 130, 160, 190 GeV 
respectively. Lines which increase toward smaller mass difference corresponds to X = a. The 
X = G lines do not show significant bino mass dependence. 

adopt simple leading order calculation, we require E 1 > 10 GeV and E T > 10 GeV. 
We define the ratio 

R(X) = r(f -> X 7 r)| aftcr cut /r(f -> Xr) (26) 

The dependence of R on the MSSM parameter is quite different between X = a and X = 
G. While f — > ar is one loop process controlled by the parameter £, the three body 
decay contains a tree level contribution which depends on non-renormalizable axino--B-gauge 
coupling (top left of Fig. Q3]) . 

When £ is small, the tree level contribution plays a dominant role in the three body decay 
into axino. This can be seen in Fig. [161 where the ratio R(a) is plotted as a function of 
the mass difference — m?- We also fix nix = 1 GeV but R(X) is insensitive to nix- 
R(a) is enhanced when the mass difference between and m f is small relative to the f 
mass. This is because the two body decay of axino is suppressed by m?~. The branching 
ratio ranges from 8% to 0.5% for the model parameters given in the figure. On the other 
hand, the branching ratio Br(f — > T'-fG) is well below 1% for the parameter given in FigfTBl 
Br(f -»■ r-fG)/Br(f -> rG) ~ 0.56% (0.84%) for m f = 100(190) GeV respectively. 
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FIG. 17: The statistics required to find the branching ratio R(a) = T(f — * vya)/T(f — > to) 
deviates more than 3 sigma from the prediction for R{G). The lines corresponds to different stau 
mass as in Figll61 

When the three body decay branching ratio turns out be above 5%, 200 stopped f is 
enough to see the 3 a deviation from the gravitino assumption. We estimate the number 
of events Advent stopped in the detector, which is required to find 3a deviation from the 
gravitino scenario as follows; 

M {BHj -> r 7 5) - Br(f -> r 7 GQ) 2 

Ntm « Br(r r 7 a) = 9 ' (2?) 

The A/event as a function of — nif is given in FigfTTl Each curve increases as — rrif 
is increased up to the value where Br(f — > r 7 a) coincides with Br(f — > tjG). When 
nig — m,f < 20 GeV, (9(1000) ( (9(10000) ) stopped f are enough to see the deviation from 
the measurement of R for £ = 0.5(1) respectively. 

Not only the branching ratio, but also the decay distribution contain the information of 
the invisible particle. The axino decay distribution can be enhanced at the region where 
x 7 = 2E 1 /m f is large and cos# is small, namely hard photon and r is back to back Q]. 
This occurs when £ is small and relative importance of the direct axino-bino-gauge boson 
coupling is enhanced. 

In Fig. [THl left panel shows the E 1 distribution for different cos# intervals. Here we 
fix£ = 1, nig = 130 GeV, and nif = 100 GeV. The distribution has significantly hard 
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FIG. 18: The distribution of E 1 in some 6 intervals for f — > t^X decay, left panel: X = a (£ 
right panel: X = G. m f = 100 GeV, and nig = 130 GeV. 
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FIG. 19: Statistics required to distinguish axino scenario from gravitino scenario, left for uIq = 
m~ a = 1 GeV and right for = rria = 30 GeV 

component for cos 9 < when compared with that of gravitino in the right panel. The 
enhancement of back-to-back events is a signature of an axino, and clear r and 7 separation 
is not required to distinguish them. 

We now estimate the the number of CNLSP decays in the stopper to see the 3 a deviation 
between r — > a^r and G^r from the decay distribution. We define A% 2 like function from 
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the differential decay width; 

AxLtK , m x , m Bi e, AE„ AB) = S, ^ ~ * ( ^ . (28) 

Here, rij(a) is the number of event in a z-th bin for f — > r7<2 when the number of stopped f 
is A^gen = 10 5 . Ai? 7 and A6 1 is the bin size in E 1 and We divide the E 1 and 6 1 into 8 bins 
and 3 bins respectively, for the range < E 1 < m? and < 9 < 7r; 

AE = rrif/8, A9 = tt/3. (29) 

The photon energy resolution of the detector is assumed AE^/E^ = 100% / \[E^J GeV. We 
apply the cut cos 9 < 0.866, E 1 > 10 GeV and E T > 10 GeV. On the other hand, fii(G) is 
the number of event in a z-th bin for f — > with = m a normalized so that the total 
number is same to that of axino three body decay. Namely, we do not use the information 
for R(x) in our fit. 

From Axd ist , we define iVdi s t(3er), the number of stopped f required to see the 3cr deviation 
between f — ► r7a and f — > as follows, 

N dist (3a)=N gen /(A X 2 di j9). (30) 

We show the iVa^o - ) as a function of - in Fig. [19j The sensitivity is significantly 
increased from the estimate using the branching ratio only. For — m? = 40 GeV, the 
deviation may be visible for (9(1000) events for £ = 0.5 (0(10000) events for £ = 1). 

Finally, we estimate sensitivity at our model points in section 3. For simplicity we assume 
f = tr. In general tr is mixed with fx,. The mixing angle is defined as 

ft = fx cos 9f + ffl sin 9 f . (31) 

For model points discussed in Table [Til the angle is sin^f ~ 0.9. The effect of the mixing 
angle in the axino decay is small because the amplitude of t l is suppressed by both by the 
small cos 2 9f factor and smaller hypercharge, and can be safely ignored. 

In FigEO the expected sensitivity at the stopper-detector is shown. Here, long-dashed 
(long-short-dashed, dashed) line corresponds to the required statistics for £ = 0.5 (0.75, 1) 
for different A, while the upper and lower solid lines correspond to the number of stopped 
CNLSP for 300 fb _1 and 3000 fb _1 . One can address the difference between axino and 
gravitino for A ~ 55(65)TeV or less for £ = 3000 fb _1 and £ = 1 (0.75). Note that the 
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FIG. 20: Solid lines are the expected number of stopped events for 300(3000)fb _1 luminosity for 
GM points with 40TeV < A < 80TeV. Dashed lines are number of required stopped events to see 
3-0- deviation in the E 1 — cos 9 distribution of f — > rja from those expected for the decay f — > tjG. 

expected reach in A is essentially determined by the SUSY production cross section, as they 
decrease steeply with increasing A. We therefore show the gluino masses scale on the top of 
the figure. 



VI. DISCUSSION 

In this paper we investigate the physics of the long-lived charged next lightest SUSY 
particle (CNLSP), which may be explored at a massive stopper-detector placed next to the 
CMS detector at the LHC. We assume the CNLSP is the lighter scalar tau, f, which decay 
into tX where X is a invisible particle. A natural candidate of the particle X is a gravitino 
but we also consider the case where X is axino a. 

In this paper, we assume very large stopper-detector next to the CMS detector, with 
total mass of 8 kton. The stopper must have a capability to measure the position where the 
CNLSP stopped, and also the energy of the r decay products. If the detector can be highly 
segmented, it is also possible to identify the r decay products separately. The number of 
stopped NLSP ranges from (9(10 4 ) to (9(100) for gluino and squark with mass below 2 TeV. 
If the size of the detector should be smaller, the number of the events must be scaled down 
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linearly. 

We estimate the statistical error for E T , which can be determined from the end point of 
the tau jet energy. We assume that the energy resolution of the detector is 150%/ ^E/GcV. 
SE T /E T ~ 60%/ 'y/ 'E T /GeV (15%/ ^E T /GeV) can be achieved for the statistics of the order 
of iVstoppcd = 100(1000). When LSP mass is above 0.2mf, one can constrain the gravitino 
mass both from above and below. From the lifetime and gravitino mass measurement, the 
supergravity Planck scale M pl can be measured. If many CNLSP can be accumulated, one 
can even study the three body decay of the CNLP f — > r^X. When the dominant decay 
mode of the f is f — > ar, we may be able to see the deviation of the decay branching ratio 
and the distribution from the X = G case. 

Finally we comment on the strategy to proceed this experiment. SUSY particles will 
be found in the early stage of the LHC experiment if SUSY scale is 0(1 TeV). If LSP 
is gravitino or axino and the NLSP is charged and long-lived, it would also be recognized 
easily. The detector proposed in this paper may be placed after the existence of long-lived 
CNLSP is observed, roughly at the same time to the high luminosity run of the LHC, or the 
proposed super LHC run. 

LHC experiment has a great potential to explore new physics in TeV regions. It is 
important to explore new possibilities that can be done with LHC. In this paper we have 
proposed a large additional detector in the CMS cavern when the long-lived CNLSP is 
found. This requires a significant modification of the CMS experiment. The reward is low 
systematics study of the CNLSP decay which primarily serve for the study of the gravitino 
sector. The determination of gravitino mass either from the lifetime or (independently) from 
E T measurement would give us a direct information of the the total SUSY breaking scale. 
Together with high precision determination of the MSSM sector expected with the CNSLP 
momentum information, the nature of the interaction of the MSSM sector and hidden sector 
can be studied in detail. We hope this paper is useful for further, and more realistic studies. 
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VII. APPENDIX 



In this appendix we list the differential decay width of tr decay into yrX where X is 
either gravitino or axino in limit where mx <C m^. Formulas in this appendix are obtained 
by taking massless limits m^/m^ — > of the formulas in Ref. Q]. The f R decay width to 
t^/G is given as follows; 



d 2 T(f R -> r 7 GQ 

dx^d COS 9 



512tt 3 [1 - (ar 7 /2)(l-cos0)] 



spins 



(32) 



where 



and 



Here 



2 _ 87r « m f rpG 



spins 



3 M*Aq 



F m (x 7 ,cos9,A B ). 



7G 

diff 



1 + cos 9 
1 — cos# 



+1 + 



Ax, 



2(1 - cos 1 



1 — x~ 



+ 



4 - Ax, 



2-x 7 (l-cos#) {2-x 7 (l -cos#)} 



1 — X, 



2 — x 7 (l — cost 



-4 + 2x - 



4(1 -x. 



{2-t 7 (1 -cosfl)} 2 



{x 7 (l + cos6») - 2A B + A B x y (l - cos6»)} 2 

4 



2A|(l-x 7 j <jx; 



2^ + 



-(x;-2x 7 -2)cosfl 



+ A B (1 + cos0)(4 — 3x 7 + x 7 cos9)(x — x, ( — 2 



— (1 + cos6 l ) 2 x 7 (x^ — s 7 — 2) 



(33) 



(34) 



Xry = 2Ey/nvf, A B = {m B /rrif ) 2 , Aq = (mg/nif)' 



(35) 



and 9 is the angle between r and 7. 

For the case of massless axino, we find 

d 2 Y{f R — > ryd) m f x 7 (l — x y ) 



where 



<ix 7 (i cos 9 



512vr 3 [1 - (x 7 /2)(l - cos 1 



spins 



^|^(^-r 7 a)| 2 



(36) 



~\|2 



spins 



ttcos 4 ^ fl 



F dm (x 



^,cos9,A B ), 



(37) 
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Xry, COS 9, Ajj) 



x: 



[1 - x 7 )(l + cos#){l + cos6» + Ag(l - cos6>)} 



+A 



{x 7 (l + cos#) - 2Az + A^x^l 
Ag(l + cos6) 

x 7 (l + cos£ 
1 . 9 1 + cos 6 



cos^)}' 



— 2v4b + AzxJl — cos( 



+^A 



cos/* 



— + — 



3a ^, f f a 

A = — ^r^°s(- 

7T COS z &w V 171 
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